Consequences of coarse grained Vlasov equations 
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The Vlasov equation is analyzed for coarse grained distributions resembling a finite width of test- 
particles as used in numerical implementations. It is shown that this coarse grained distribution 
obeys a kinetic equation similar to the Vlasov equation, but with additional terms. These terms 
give rise to entropy production indicating dissipative features due to a nonlinear mode coupling 
The interchange of coarse graining and dynamical evolution is discussed with the help of an exactly 
solvable model for the selfconsistent Vlasov equation and practical consequences are worked out. 
By calculating analytically the stationary solution of a general Vlasov equation we can show that 
a sum of modified Boltzmann-like distributions is approached dependent on the initial distribution. 
This behavior is independent of degeneracy and only controlled by the width of test-particles. The 
condition for approaching a stationary solution is derived and it is found that the coarse graining 
energy given by the momentum width of test particles should be smaller than a quarter of the kinetic 
energy. Observable consequences of this coarse graining are: (i) spatial correlations in observables, 
(ii) too large radii of clusters or nuclei in self-consistent Thomas-Fermi treatments, (iii) a 
structure term in the response function resembling vertex correction correlations or internal structure 
effects and (iv) a modified centroid energy and higher damping width of collective modes. 



I. INTRODUCTION 

The Boltzmann kinetic equation is a successfully used 
dynamical model to describe heavy ion collisions and 
multifragmentation up to intermediate energies [[j] pd| . 
The scattering is added to the Vlasov equation either 
by the relaxation time P4-|lq] or by scattering integrals 
of the Boltzmann-type p|,[7[|16|-|2^] or their nonlocal ex- 
tensions p3|-p5t in the same way as it was done in the 
Landau- Vlasov equation for hot plasmas. All these sim- 
ulation have in common that the motion of particles be- 
tween collisions are covered by Hamilton equations or the 
Vlasov kinetic equation with a selfconsistent potential. 

Even within collisionlcss Vlasov codes the onset of frag- 
mentation is described p6|~p8[- The use of selfconsis- 
tent Vlasov equation is not restricted to nuclear colli- 
sions but has been applied successfully for collisions of 
ions with metal clusters [p9| Despite the fact that 
the Vlasov equation is a reversible kinetic equation and 
the initial configuration should be retained after a long 
enough calculation, fragmentation and energy dissipation 
is observed. Obviously this is due to the fact that the 
Poincare time is much larger than the time where the 
phase space is filled by various trajectories. Therefore 
one can consider superficially this spreading as an irre- 
versible process in the sense of entropy production. Nev- 
ertheless, the underlying dynamics is reversible. 

The numerical implementation of Vlasov codes de- 
mands a certain coarse graining of the space and momen- 
tum coordinates. This numerical uncertainty is quite suf- 
ficient to generate genuine dissipation and entropy pro- 
duction. This fact has been investigated in J32| and will 
be considered in detail within this paper. It has been ar- 



gued that the errors due to numerically coarse graining 
accumulate diffusively. We will demonstrate that these 
errors will lead to a unique equilibration of Vlasov equa- 
tion. While the diffusion itself is less effected, the coarse 
graining leads to a different dynamical evolution and has 
consequences for the extraction of damping rates from 
Vlasov-type simulations. 

Theoretically one can derive kinetic equations by 
phase-space averaging of the Liouville equation which 
itself is exactly of Vlasov type |3^]. The resulting col- 
lision integrals represent two different approximations 
of the nonequilibrium dynamics: (i) The truncation of 
coupling to higher order correlations (hierarchy) and (ii) 
The smoothing procedure which translates the fluctuat- 
ing stochastic equation into a kinetic equation for the 
smoothed distribution function. This latter procedure 
is sometimes also called coarse graining. For the discus- 
sion of appropriate collision integrals also in the quantum 
case see j34j] . Here we will focus on a detailed analysis of 
consequences of coarse graining. 

The idea used here dates back to the work of Gibbs 
and Ehrenfest |35||36| . They suggested to coarse grain the 
entropy definition by a more rough distribution function 



f(p,r,t) 
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A(p,r) 



fs(p',r',t)dp'dr'. 
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The physical meaning consists in the fact that any ob- 
servable is a mean value of an averaging about a cer- 
tain area in phase space. It was shown that the entropy 
with this coarse grained distribution increases fl37| . This 
means that in a closed system the entropy can rise if we 
average the observation about small phase space cells. 
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The interpretation is that other phase space points can 
enter and leave the cell which is not compensated. A 
phase space mixing occurs J37[ since the two limits cannot 
be interchanged, i.e. the thermodynamic limit and the 
limit of vanishing phase space cell. The coarse graining of 
Ehrcnfest can be observed if the thermodynamical limit 
is carried out first and the limit of small phase space 
cells afterwards. It solves therefore not the problem of 
entropy production, but gives an interesting aspect to 
entropy production by coarse grained observations [p7fl . 

In this paper we like to investigate three questions: 

• Which kinetic equation is really solved numerically 
if the Vlasov equation is implemented in numerics? 

• What are the properties of this kinetic equation, 
especially which kind of dissipative features appear 

? 

• What are the consequences to practical applica- 
tions, e.g. damping of giant resonances and binding 
energies? 

The outline of the paper is as follows. Next we de- 
rive the kinetic equation which is obeyed by the coarse 
grained distribution function. Then in chapter III we dis- 
cuss the entropy production. We demonstrate with the 
help of two exactly solvable models that this entropy pro- 
duction is due to mixing, i.e. a mode coupling and not 
simply by spreading of Gaussians. The solution of the 
stationary Vlasov equation is then presented in chapter 
IV. We will find that the stationary solution can be rep- 
resented as an infinite sum of modified Boltzmann dis- 
tributions. This expansion shows the unique character 
of time evolution which is only determined by the ini- 
tial distribution. In chapter V we discuss consequences 
of this result: (i) The thermodynamics becomes modi- 
fied by spatial correlations. The selfconsistency leads to 
a modified Thomas- Fermi equation lowering the bind- 
ing energy, (ii) the structure factor shows a substructure 
similar as obtained from vertex corrections and (iii) the 
damping width of collective resonances is shown to be 
larger by coarse graining. While the centroid energy is 
smaller by momentum coarse graining it increases by spa- 
tial coarse graining. 



II. COARSE GRAINED VLASOV EQUATION 

The origin of the coarse graining may be the numer- 
ical implementation or the use of averaged distribution 
functions instead of the fluctuating one. To illustrate 
the method we examine the quasi-classical Vlasov equa- 
tion and show which equation is really solved if one is 
forced, by numerical reasons, to use coarse graining. It 
will become clear shortly that instead of Vlasov equa- 
tion a modified kinetic equation is solved when coarse 
graining is present. The quantum mechanical or TDHF 



equation can be treated in analogy. We start from the 
Vlasov equation 

^-Jsiprt) + --^-fsiprt) - ^-V s (r,t)^-f s (pH) = 
ot mar Or op 

(2) 

which solution can be formally represented as an infinite 
sum of exact test-particles 



f s (prt) = ^(*)) S(P - m)) 



(3) 



where the test-particle positions and momenta evolve cor- 
responding to the Hamilton equations Ri (t) = Pi/m and 
Pi{t) — —dRV(Ri(t)). In the following we understand 
p, r, P and R as vectors and suppress their explicit nota- 
tion. 

For the mean-field term V$ we assume a Hartree ap- 
proximation given by a convolution of the density with 
the two-particle interaction V 



V s (r,t) = / dr'V (r-r') 



dp' 



-fsip'r't). (4) 



In practice all numerical calculations use two assump- 
tions : (i) The infinite number of test-particles is trun- 
cated by a finite value, (ii) The used test particle has 
a finite width due to numerical errors and/or smoothing 
demand of the procedure. While in [h| was discussed 
that the approximation (i) leads to a Boltzmann-like col- 
lision integral, the approximation (ii) will deserve further 
investigations. Especially, we will show that the finite 
width of test-particles leads to a coarse graining and a 
dissipation forcing the system to a Boltzmann- like dis- 
tribution. This is even valid with infinite numbers of 
test-particles. Therefore we consider the effect of coarse 
graining as the most determining one for one-body dissi- 
pation. 

The finite width of test particles can be reproduced 
most conveniently by a convolution of the exact solution 
(|) with a Gaussian g a (x) = (2Tral)- 1 - 5 exp{-x 2 /{2a 2 a )) 
resulting in the coarse grained distribution function / 



f(prt) = {fs} 5 



dp' 

—^dr'g r (r - r')g p (p - p') fs{p'r't) 



(2-Kh) 

OO 

= J29r(r-Ri(t))g p (p-Pi(t)) 



(5) 



To answer the question which kind of kinetic equation 
describes this smoothed distribution function, the kinetic 
equation for ^ is derived from (||) by convolution with a 
Gaussian. The equation for general coarse grained mean- 
fields has been already derived in (io). In order to make 
the physical content more explicit we calculate the differ- 
ent terms explicitly and present the coarse grained kinetic 
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equation. The free drift term ^d r fs takes the form after 
convolution 

{-d r fs} a = ^f(prt) + —4-4-fiprt) (6) 
to mor to or op 

which is established by partial integration. We see that 
the free streaming is modified by an additional resistive 
term which will give dissipative features. 
The mean-field term takes the form 



i) 







d r d 



{gZ V s(r)—fs(prt)} g = —{ — V s (r)f s (prt)} 



dp 



dp di 

= ^{{§- r Vs(r)fs(prt)} 9r } gp . (7) 



The space convolution with g r is performed using the 
relation Q 



{AB} g = {A} 9 expK 2 d r d r ){B} g 
with the result 



(8) 



g-{g^V(r,t)exp(a 2 r d r dr){fs(prt)} gr } gp . (9) 



V(r, t) is the mean-field calculated with the space and 
momentum coarse grained distribution / instead of fs, 
via (g) which can be seen as 



V(r,t) = {V s (r,t)} gr = J dr"V (r") J dr' g r {r - r" - r') 
dp' , ( , , ,s 

dr"V (r") J dr' g r {r - r" - r') 

dp' f dp , , , 

9 P \P~P )J8(P ,r ,t) 



{2-Khf J {2irfr) 3 
dr"V (r") ' ''>' 



:f(p,r-r",t) (10) 



where the last equality shows the invariance of particle 
density due to coarse graining. The momentum convolu- 
tion with the Gaussian is then performed in (^|) to yield 
the momentum and space coarse distribution function /. 
We like to point out that the test-particle method would 
lead to a further folding of the mean-field potential if it 
is read off from a finite grid [Q . 

The coarse grained Vlasov equation reads now 

ot to or to or Op 

Q ^ ^ Q 

- ^-V{r,t)exp(a 2 d r dr)^-f(prt) = 0- (11) 
o r op 

This equation is the main result of this chapter and repre- 
sents the Vlasov equation for the coarse grained distribu- 
tions and should be compared with the Husimi represen- 
tation of GfJ. While the distribution function is exactly 



the Husimi representation of the Wigner function the cor- 
responding coarse grained kinetic equation was not given 
before. 

Equation (^) represents the kinetic equation which 
is really solved numerically if the Vlasov equation is at- 
tempted to be solved. The coarse graining leads to two 
additional contributions besides the original Vlasov equa- 
tion. We will see that this causes just the dissipative like 
features. While we will continue now and investigate ( fTl| ) 
more closely, in appendix [a| the question is addressed how 
the underlying dynamics is modified by coarse graining. 
We find that in principle the coarse grained equation can 
be mapped to the original Vlasov equation if one defines 
new testparticles with modified dynamical equations by 
a modified potential. The relation between the original 
and the modified potential appears just as the inverse 
folded mean field potential which is the opposite relation 
than presented in 



III. ENTROPY PRODUCTION 



To analyze the dissipative feature we rewrite ( |11| ) into 

i-J(prt) + 2~%-f{prt) - ^-V(r,t)SLf(prt) = I diss 
ot mor or op 



(12) 

with the one-body (collision integral like) term 

T ^ 9 9 ft A 

mor op 

Q ^ ^ Q 

+ ^-V(r,t)(eM^rdrdr) - l)^-f(prt) 
Or Op 

a (^V(r,t) - |) A A /(prt) + o(a 4 ). 
or z m or op 

(13) 

In order to demonstrate the entropy production ex- 
plicitly we use the linearized form of (113) and built up 
the balance equation for entropy S — /In/ by multiply- 
ing (^) with In/ and integrating over p. The entropy 
balance reads then 

a, ,n . 9 f dp p f dp d r fd p f 

(14) 



dr J (2-kK) 3 m" 



with $(r,t) = -cr^^2-V(r, t) + The entropy density 
obeys therefore a conservation law with a source term on 
the right hand side of (|l4|). Especially one sees that the 
total entropy change is equal to 



S(t) 



drS{r, t)= dr$(r, t) 



dp d r fd p ! 



(2nh)3 f 



(15) 
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This expression allows already to learn some interesting 
properties of such distribution functions which lead to 
an entropy increase. For a distribution function symmet- 
ric either in p or r, no entropy production occurs since 
df would be asymmetric. We get an entropy produc- 
tion only for explicit space dependent distributions due 
to the d r f factor. Near equilibrium where the assumed 
nonequilibrium distribution functions are falling in space 
and momentum, the second derivative of the mean-field 
is negative (due to finite systems) and therefore $ > 0. 
Consequently we obtain an increase of entropy on average 
due to the spatial coarse graining. This means that the 
finite width of test-particles induces an entropy increase 
similar like irreversibility. 

The reason for the entropy increase is a nonlinear mode 
coupling which can be considered as a general feature of 
coarse graining. Therefore we rewrite the collision inte- 
gral in the first line of ( |l3| ) by Fourier transformation into 
another form 



J dr" VK 3 J (2tt 



S(r-r')fr d r -^-V(r,t)) 
m or 

dk\dk 2 



nf 

x exp(— i(ki + k 2 )r — fj^fc 1 fc 2 + ik 2 r' + ik\r")\ . 

(16) 

We see that due to the spatial coarse graining a nonlinear 
mode coupling occurs. This is represented by the product 
kik 2 between the modes in the distribution function and 
the mean field. This latter effect is the reason for the 
production of entropy and is connected with the spatial 
coarse graining oy.[] 

To understand the physical origin of this entropy pro- 
duction we choose two simple models for illustration. In 
the first example we assume a fixed external potential. 
In the second example we give an exactly solvable model 
including the selfconsistent mean-field potential. 



A. Free Streaming 

The initial condition before folding is assumed as 

fs(rp0) cx e'^5{r). (17) 

The Vlasov-equation with V — yields then the time 
dependent solution 



* Please remark that / and V are coarse grained values it- 
self which Fourier-transforms contain an exponential kia r /2 
and fc|a^/2 respectively. Therefore the total exponent (fci + 
&2) 2 0"r/2 appears and the expression is convergent while the 
bare fcifc2 product above may lead to the impression of non- 
convergence. 



fs(r,p,t) <x e~*> 6(r-^t) 
m 

or in vector notation x = (r, p) 



fs(r,p,t) 



2tt 



3/2 



(£) 



3/2 



e -A^_ f(r _ pi) 2 



ji — ► oo 



(18) 



(19) 



(2 7 r) 3 v /Dct{A} 



with 



1 -fit X + fit 2 



The distribution ( [T9| ) is to be folded with 
G{r,p) = 



(20) 



e - 5 x i, J 



27r) 3 A/Det{S} 



C7~ A 

o 



Using the Gaussian folding theorem, we obtain 



f{r,p,t) 



e 2 T 



(21) 



(22) 



(2ir) 3 \/ Det{M} 
M = A + S 




and the entropy is 

S=- J drd P flog(f)= J drdpf(r,p,t) 

-xM _1 x + 3 log (2tt) + - log ^Det{M}) 



(23) 

Using known Gaussian integration rules and performing 
the limit /i — ► oo, we obtain finally 



^ = 3(log(2 7 r) + l) + flog(^ + ^ f + ^ 



(24) 



which appreciates S(t) = 3 log ^ ^J e tj for large times. 

Interestingly the long time limit is only modified by a p . 
One can see that the entropy is monotonically increas- 
ing. The reason is the continuous solvation of the folded 
distribution. 



B. Selfconsistent bounded model 

After demonstrating the increase of entropy with time 
we like to consider two questions: (i) Is this increase due 
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to the spreading of the Gaussian, which we assumed for 
space coarse graining ? (ii) Can we interchange the pro- 
cedures coarse graining and dynamical evolution? This 
means we like to check whether we will get identical re- 
sults when we first solve the exact Vlasov equation and 
then coarse grain the solution or when we coarse grain 
the Vlasov equation first and solve the modified one af- 
terwards. This question reveals the sensible dependence 
on the initial distribution. 

To answer both questions we consider another example 
of exactly solvable model given in [Q , where we will add 
an external harmonic oscillator potential. The poten- 
tial V(r, t) consists of the separable multipole-multipolc 
force 1*1234 = w<?i2<?34 for the mean-field and an external 
harmonic oscillator potential 



1 



2„2 



V(r,t) = vg(r)Q(t) + -muj z r 

where the selfconsistent requirement is 
drdp 



Q(t) 



T9(r)fs(r,p,t) 



(25) 



(26) 



according to (||). We may think of the external har- 
monic potential as a representation of a realistic con- 
fining potential expanded around the origin U cx t{f) ss 
£4xt(0)(l - (r/i?) 2 /2) with the radius R. We have there- 
fore 



« -{/ cxt (0)/mi? 2 



(27) 



The one-particle distribution function obeys the quasi- 
classical Vlasov equation (^|). We can solve the Vlasov 
equation exactly by solving the differential equation for 
the equipotential lines, which are the Hamilton equa- 
tions. We choose a form factor of 



g(r) = a x r x + a y r y + a z r z 



(28) 



This model is special by two reasons. Firstly, for this 
model the linearization (|l^) is exact, because higher than 
second order spatial derivatives vanish. Secondly, within 
this model the quantum- Vlasov equation agrees with the 
semiclassical Vlasov equation investigated here. 

The Hamilton equations of trajectories which corre- 
spond to this model 



are solved as 



dtp = —vaQ(t) — mu) 2 r 

a P 
o t r = — 

m 



sin Lot cos tot 

ma) cos Lot —mco snoot 
t 

dt'Q^f^^-p 

K ' \ COS LO(t — t ) 




(29) 



(30) 



Now we know that the constants of motion c\ , ci are con- 
stant at any time of evolution. Therefore we can relate 
them to the initial momenta and positions which gives 



Po 

ro 



(31) 



From equation ( |30|) we express now ro and po as functions 
of (r, p, t) which results into 



cos Lot — — sin Lot 



mLO sin Lot 



cos Lot 







(32) 



Given an initial distribution /o(ro,Po)i we can express 
the solution at any time as 



fs(r,P,t) = /oOo,Po) = /o(Ax + C) 



(33) 



with x = (r, p) and the matrix A and the vector C can 
be read off from (|32"|). The selfconsistency requirement 
(p6|) leads to a determination of Q(t) 



Q(t) — a < r > cosh\/ Q 2 — u) 2 t 



;\ZQ 2 - Lo 2 m 



< — > sinhv^ 2 - lo 2 t (34) 



where 



n 2 

N 

<r > 



<^> 
m 



ro 



-N 

m 

drdp 



(2ttH) 3 
drdp 



(2-Kh) 
drdp p 



fo(r,p) 
3 rf {r,p) 



fo{r,p) 



(35) 



are expressed by moments of the initial distribution. 

We have given an exact solution of the selfconsistent 
3-D Vlasov equation. It is interesting that for £1 > lo 
we have a positive Lyapunov exponent \/£l 2 — io 2 while 
in the opposite case the solution oscillates. We like to 
point out that ([33]) represents a general solution of any 
Vlasov equation if we understand Ax + C as a nonlin- 
ear transformation solving the corresponding Hamilton 
equations. 

In analogy to the foregoing chapter we can now coarse 
grain this exact solution. We are able to check by this 
way whether the so far observed entropy increase is due to 
the spreading of the test particles with Gaussian width. If 
so, we would observe an entropy increase even if we start 
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with an equilibrium solution. We will demonstrate now 
that this is not the case. Instead, we have a nonlinear 
mode coupling. The actual calculation is performed in 
the appendix pi We choose as an initial distribution a 



Gaussian one. The entropy finally reads from (B7) 



,9 = 3 log (27re) + log • m'Kc 2 ^ ) i 

+ {^- + <Jp<Jr) + \{m 2 u 2 a 2 r -<7l){ x " --)cos{2u>t) 



1 



^_ 1 

Xm 2 u 2 (i 
1 



(36) 



The entropy is just oscillating around a stationary value. 
The Gaussian form of initial distributions leads to no 
stationary solution. This shows that the entropy increase 
observed so far is not due to the spreading of Gaussians 
but due to the form of initial distribution. Other initial 
distributions would lead to an increase of entropy which 
can be checked e.g. with ground state Fermi functions. 

Now we turn to the second question concerning the in- 
terchange between coarse graining and solution of kinetic 
eguation. While we have solved the exact Vlasov equa- 
tion so far and coarse grained the solution afterwards we 
now revert the procedure and solve the coarse grained 
Vlasov equation ([l2|) with (|l3|). Please remind that the 
linearization with respect to the width is exact for this 
model here. The entropy in this case is calculated from 



( |BllD with the result 



5 = 31og2 7 re + ^log{i- + i(a 



2 2 2 2\ 

p -TO W 0~ r ) 



x ( 



1 



\m 2 u) 2 



1 
/' 



oy)(l — cos(2a;t)) 



(37) 



Comparing with (|36|) we see a different expression, how- 
ever the oscillating behavior remains. Within linear 
orders of a both expressions differ by a constant of 

2-p- + 2^. The difference between these two expressions, 
which corresponds to the interchange of coarse grain- 
ing and dynamics is explained by the use of the same 
initial distribution in both cases. We obtain a differ- 
ent dynamical behavior. If w e wou ld additionally coarse 
grain the initial distribution ( B10| ) we would have to re- 
place AbyA + S" 1 . Instead of ( |B12| ) we would obtain 
M = (A T (A + S~ 1 )A)~ 1 — B and the resulting entropy 
agrees exactly with (pq). As already pointed out this 
would be not the case for nonlinear models other than 
( p8| ) since then the selfconsistent condition (^6|) will be 
affected by coarse graining itself. 



C. Consequences 

This observation has some practical consequences. 
Since in all numerical calculations one solves the coarse 
grained Vlasov equation instead of the exact one, which 



corresponds to first coarse graining and then solving, one 
should not expect the correct dynamical behavior start- 
ing from a fixed initial distribution. Instead the correct 
procedure is to coarse grain first the initial distribution 
with the Gaussians of the used test-particles and then 
solve numerical the time evolution of this distribution. 
The refolding at any time step yields then the exact so- 
lution of Vlasov equation for such trivial models as de- 
scribed here where the selfconsistency condition is not 
altered by the coarse graining. 

It is important to stress again that the above equiva- 
lency between the two ways of solving and coarse grain- 
ing are only equivalent in the linear model ( p8[ ) which 
we have solved exactly. For more realistic potentials the 
selfconsistency condition like (|2^) becomes altered itself 
by the coarse graining which leads to an essential non- 
trivial change in the dynamics. Therefore both method 
will lead to different results. The nonlinear dynamics is 
non-trivially changed if one coarse-grains first and solves 
then or if one solves first and coarse-grains afterwards. 
However, in order to diminish the numerical error due 
to coarse graining the refolding should be performed at 
least. 



1. Time scale of entropy production 

From the chapter above we see that the entropy pro- 
duction due to coarse graining in an harmonic external 
potential is oscillating like 1 — cos (2ojt). We can use this 
fact to define the typical time scale of entropy production 



2lj 



t/ext(0) 



(38) 



where we used (|27|). For a typical nuclear situation we 
obtain therefore r c w 4.8A 1 ^ 3 fm/c such that for 12 we 
get t c w 12fm/c and for 2Q8 Pb we have r c w 28fm/c. It 
is remarkable that the time scale on which the entropy 
is changing is independent of the used internal interac- 
tion. This clearly underlines the Landau damping type 
of dissipation. Coarse graining produces no genuine dis- 
sipation. 



IV. STATIONARY SOLUTION OF COARSE 
GRAINED VLASOV EQUATION 

We will demonstrate in this chapter that the solution of 
the general coarse grained Vlasov equation is approach- 
ing a modified Boltzmann limit for long times provided 
a stationary solution is approached at all. The example 
in chapter |IIIB has shown that this is not the case ev- 
ery time. Exclusively for appropriate initial conditions, 
which will be characterized in chapter IV C , we will ob- 



tain stationary solutions due to the additional terms in 
the Vlasov equation (jlj) . 
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A. Solution of coarse grained Vlasov equation 

We want to consider the solution of the kinetic equa- 
tion ( |l2| ) with the assumption of arbitrary time depen- 
dent mean-fields. Therefore we solve the following partial 
differential equation neglecting the selfconsistcncy in V 
at first and remember it later. We consider 

^r/(prt) + *(r, t)£~£-f(prt) + -^-fiprt) 
at drop mor 

-^ yM) ^ /(prt)= ° (39) 

with <f>(r, t) = (-^■ — a^S^V(r,t)) and the boundary con- 
ditions /(r,p,0) = f (rp),J dpdrf(rpt)/(2irh) 3 = N. 

This equation is of parabolic type and as an initial 
value problem well defined with an unique solution. The 
unique stationary solution, if such a stationary solution 
is approached, is an implicit representation of the sta- 
tionary solution by the remembrance of the dependence 
of V on the distribution function itself, i.e. the selfcon- 
sistency. Because the problem is uniquely defined it is 
enough to find a special representation of the solution. 



1. Representation of stationary solution 



The stationary solution of (|3£|) can be found by sepa- 
ration of variables. Assuming / s tat(p^) = fp(p).fR( r ) we 
have 



iSatbO = const x exp 



-c (— 
2m 



, d'V(r') 
dr' r K ,J 

(40) 



with the separation constant c n which holds for vectors 
p, r. The general stationary solution is given by superpo- 
sition of these c„ dependent expression ( fief ) . Linearizing 
via the coarse graining a leads us to 



/stat(p,r-) = y^a»< 



-Cn(fer+V(r)) 



.<J 2 __. . a 2 



x (l- c 2 n (^V(r)-^(V'(r)) 2 )) (41) 

from which we deduce the equilibrium density distribu- 
tion to have the form 



2Trh 2 c r 



)3/2 e - C „V(r) 



2 2 

x { l- c 2 n (^V(r)-^-(V'(r)) 2 )), (42) 
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Please remark that the summation can also be replaced 
by an integration which would translate into continuous 
functions c„ and a n . For the reason of legibility we dis- 
cuss further on only the discrete case. 



With this solution we have derived our goal to present 
a stationary solution of the modified Vlasov equation. 
Because the initial problem was uniquely defined this so- 
lution is the unique stationary one. It has to be remarked 
that due to the selfconsistency ( dependence of V on the 
distribution function itself ) equation (|4l| ) is an implicit 
representation of the solution. 



B. Determination of c„ 

The open expansion coefficients c„ are completely de- 
termined by the initial distributions. This can be seen 
as follows. To solve the Vlasov equation we can rewrite 
the solution into Hamilton equations with a time depen- 
dent potential. This time dependence comes from the 
selfconsistent potential and/or from an external poten- 
tial. The selfconsistency is then represented by a nonlin- 
ear determining equation similar to (|2^). The Hamilton 
equation for the trajectories can be solved in principle 
with integration constants Cj as we demonstrated in the 
model. Because these are constants at any time we obtain 
a transformation between initial coordinates and the co- 
ordinates at a later time, including the time dependence 
as (r ,p ) = ^Kp, £]• Since we like to solve the Vlasov 
equation as initial value problem with a given initial dis- 
tribution fo(ro,po), the time dependent solution is given 
formally by f(r,p,t) = f (A[r,p,i\). //this solution ap- 
proaches a stationary state, which is dependent on the 
initial distribution as well as on the interaction poten- 
tial, we obtain with the result of the foregoing chapter 



fo{A[r,p,oc]) = ^2a n e~ 



2 2 

x (i-4(^F(r)-^(F'(r)) 2 )). (43) 

777 2 

By integrating over r and inverse Laplace transform it 
is possible to extract the c„ uniquely. In the next chap- 
ter we will find that only such initial distributions will 
lead to stationary solutions which obeys c„ < ^ • Since 

the c n are determined by the initial distribution and the 
used potential, we can decide which initial distributions 
will lead to stationary solutions with a given interaction 
potential and coarse graining a p . 



C. Global stability of the stationary solution 

Since the general possible solution of the stationary 
coarse grained Vlasov equation covers a range of unphys- 
ical solutions which will never be reached, we have to 
ask which solutions are the stable ones. Therefore we 
employ a linear response analysis of the Vlasov equation 
( |39| ) to an external potential U ex t- We assume a ho- 
mogeneous equilibrium characterized by the distribution 
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/o (p) ■ Then the equation ([59]) can be linearized according 
to f(prt) = f (p) + 5f(prt) as 

(—iLU + i — )5f(pkui) — iV' (n )kdr,fQSn(k,uj) 
m 

+ ikdpf U ex (ku;) 

- i^-kd p df( P kuj). (44) 

Here we have Fourier transformed r, t coordinates. It can 
be easily solved for Sf. Integrating this solution over p 
an algebraic solution for Sn is obtained 



Sn = U„ 



n(few) 



1 - V {n )Tl(kLu 
with the polarization function 
dp 



(45) 



n(few) 



{2nh) 3 
dp 



(1 



p kdp kdpfo 



pk 



m oj 

(2tt^) 3 ^ ' m 2 (cc-^) 2 



pk 



(1 



1 



kdpfo 



(46) 



We see that the usual RPA response function is modified 
by a structure function 



M 2 {kpuo) = 1 + 



fcV 2 



1 



(to- 



pk\ 2 ' 

m. > 



(47) 



This describes the fact that the elementary particle con- 
sidered here (testparticles) have a finite width or an in- 
ternal structure. For large distances k — > we see that 
M approaches 1, which means that this structure is only 
resolvable at smaller distances. 

A known approach to find the structure functions M 
inside the RPA polarization function is to include ver- 
tex corrections. The resulting response functions can be 
written generally 



dp M 2 (kpu) Wpfo 



{2nh) 



(48) 



where d n is the spatial integral over the potential depen- 
dent part. From ( f46|) we obtain as criterion of stability 



„ , k 2 G 2 d 2 , m 2 uj x - 
Imn = - 1 + —2.—)-— 

m z ouj z 2-kk ' 



> 0. 



(50) 



To provide stable solutions we demand that any term in 
the sum should be positive. We observe now two sources 
of possible instability: (i) The expansion coefficient a n 
without coarse graining can become negative for enough 
pathologic initial functions. For such initial distributions 
we would not get stable solutions by Vlasov dynamics 
itself, (ii) Provided the original distribution is stable 
a n > we find an additional criterion for stability if 
we use coarse graining. Demanding the coefficients to be 
positive leads us from (f4rj|) to the equation 



1 



2 / 3c n 
°p( 



•<(|r)>o. 



(51) 



In the case of small wave lengths 
most restrictive condition [j] 



k 



< 



2a 2 p 



oo we obtain the 



(54) 



If this condition is not fulfilled we find muj/k combina- 
tions such that the sum in (^o|) can change the sign and 
the solution is unstable. In other words ( |54| ) is a nec- 
essary condition for stability of the stationary solution. 
Since the c„ are determined by the initial distribution 
acc ordin g to (|4^) and the potential as discussed in chap- 
ter IV B , we see that only special initial distributions can 



lead to stationary long time distributions due to coarse 
graining with a given potential. 

Interestingly, for a Maxwellian initial distribution c„ = 
1/fcgT and a n — const^o, the coarse graining has to 
obey 



By approximating this structure function M by (|47|) we 
might simulate higher order correlations by finite momen- 
tum widths of testparticles. Also the influence of finite 
size can be compared with this expression 

The equilibrium solution is stable as long as the in- 
crement of the collective mode does not change the sign. 
Otherwise the collective mode would exponentially in- 
crease with time and the solution would be unstable 
with respect to small perturbations. The collective mode 
is given in linear response b y t he complex solution of 
1 — y' (?io)n(fca->) = of eq. (|4q). For small increments 
the sign of the complex part is completely determined 
by the sign of Imll of ([46]). Let us check this stability 
condition for the stationary solution. From (J40j) we have 



/stat(p) = / drf stat (p,r) = ^a„e" 



Li 7) 



(49) 



^The less restrictive condition for stability reads 



3fc 2 
2muj 2 



1 - 



4m 2 uj 2 



In the general case with arbitrary wave length we see that 
© is fulfilled if 



2 m\uj\ 
< 3 k 



(53) 



We conclude that all solutions ( |40[ ) are stable if the coarse 
graining is smaller than the typical scales in the system ([si]). 
If ( |53| ) is not fulfilled, than only such solutions (|Icj) are stable 
solutions which expansion coefficients c n are smaller than the 
inverse coarse graining via (|5^). Therefore (^) is the most 
restrictive condition. 
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Op < 



m 



(55) 



to provide stable solutions. This gives the intuitive clear 
interpretation that the coarse graining energy a 2 /2m 
should be smaller than a quarter of the kinetic energy 
in order to render numerical investigations stable. 



V. THERMODYNAMIC CONSEQUENCES 

We have reached our goal to show that phase space 
coarse graining leads to dissipative features of reversible 
kinetic equations and the special coarse graining with 
Gaussian width forces the system to a sum of modified 
Boltzmann - like distributions. We have shown alge- 
braically the numerical observed fact that in ordinary 
Boltzmann codes two different limiting values of the dis- 
tribution functions are appreciated |p9| . From one- par- 
ticle dynamics a Boltzmann- like distribution and from 
quantum Boltzmann collision integrals a Fermi distribu- 
tion evolves. In contrast to earlier work we find that 
even for an infinite number of test-particles the Boltz- 
mann limit is appreciated due to the finite width of test- 
particles. 

From the found stationary solution of the coarse 
grained Vlasov equation we can now derive thermo- 
dynamic consequences. We assume for simplicity a 
Maxwellian initial distribution which may evolve in time 
via the coarse grained Vlasov equation. Then the coef- 
ficients c n of the stationary solution ( pfl| ) are only single 
parameters determined by c„ = (3 = l/fc^T and 



= 5 n ,t 
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J dre~ 



■0V(r 



[l-P[V{r)%-°-f (d r V(rW]} 



(56) 

by the requirement of particle conservation that the in- 
tegral over momentum and space of ( ftl| ) should equal to 
the particle number N. Here g describes the spin-isospin 
degeneracy and \\ — 2ttH 2 /rnksT and we use the lin- 
earization in a since the original equation (|3^) is valid 
only up to orders of o(a ). 

From the expression (|4l]) we see that the coarse grain- 
ing leads to a modification factor of the distribution func- 
tion / s tat (p r ) m comparison with the Maxwell-Boltzmann 
distribution function fjaip^) = const exp(— f3p 2 /2m — 
(3V(r)) 

Wpr) _ l-P 2 (v(r4-4(drV(r)?) 
= i _ f rl [V (r)- < V(r) >} 



> 



\ [{d r V{r)f- < (d r V(r)) 2 >] 



where the spatial average is abbreviated as 
< a >= fdre-WWa/fdre-WW. 

As a consequence, no contribution of coarse graining 
occurs for mean values of momentum dependent observ- 
ablcs within the lowest order of coarse gaining width. 
However, for space dependent observables one gets a 
modified thermodynamics by a p coupling the mean spa- 
tial fluctuation of the potential and by a r coupling the 
fluctuation of the potential derivative. Consequently, 
only the spatial thermodynamical quantities will be in- 
fluenced by the underlying coarse graining. We can con- 
sider this modification factor as the expression of spatial 
correlations induced by momentum and spatial coarse 
graining. 
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FIG. 1. The influence of different coarse graining parame- 
ters on the selfconsistent potential. Above the space coarse 
graining is zero and below the momentum coarse graining is 
set to zero. While the momentum coarse graining leads to 
an increase of the selfconsistent potential the spatial coarse 
graining enhances the gradient and produces a skin. 



A. Selfconsistency 

Now we return to the question of selfconsistency. 
We have so far assumed silently that (f|) or ( |Io[ ) are 
completed by a selfconsistent potential V(r). With- 
out coarse graining this would be accomplished by solv- 
ing the Thomas-Fermi-like equation (^) with / cx 
exp-f3{p 2 /2m + V(r)) 



V(r) 



N J dr'V a (r - rOe-^M 



.9 



/ dr'e~ 



■PV(r>) 



(58) 



(57) 



and analogously for the degenerate case. For the coarse 
grained case we see now from (41) and (|56|) that the 
potential in the distribution function / has to be replaced 
by 
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V cS (r) =V(r)+p 



JL V (r)-f(d r V(r)y< 



(59) 



Therefore one can solve the Thomas Fermi like equation 
([58|) as before but replace afterwards V — ► V c g. In the 
following let us assume we have solved the Thomas-Fermi 
equation (^8|) and we will discuss now the influence of 
coarse graining. 

We see that the coarse graining in momentum space 
increases the potential ([59]) globally. In contrast, the 
coarse graining in space causes a widening of the poten- 
tial. This is demonstrated in figure ^, where we assume 
a Wood-Saxon potential for V(r) in 208 Pb and plotted 
the change of the potential for different parameters of 
coarse graining at a temperature of 10 MeV. While the 
momentum coarse graining leads to an overall increase, 
the spatial graining sharpens the gradient in the potential 
and causes the appearance of a skin. 

One can consider this as a modification of the Thomas- 
Fermi equation by the finite width of phase space grain- 
ing. The coarse graining obviously produces higher bind- 
ing properties and larger rms- radii. The appearance of 
a skin is important to note as a relict of the testparticle 
width. 



B. Consequences on collective modes 

As a practical application we show now that the coarse 
graining which is numerically unavoidable leads to false 
predictions concerning the energy and the width of giant 
resonances. For illustrative purpose we restrict to the low 
temperature T = case. The collective mode is given by 
the complex zeros to — E — iT of the denominator of ( 45 ) . 



This solution provides us with the centroid energy E and 
damping T of the collective mode. We find for T = the 
dispersion relation 



l-VIL(q,u) = 1 + F (1 + 



2pj du 2 ' 











X 


-H 















i-u&(l 
2 v 



(60) 



where we have introduced u = muj/pt/k, the Landau 

parameter F = 2V | r "°^ mp/ 
p f . 



and the Fermi momentum 




FIG. 2. The centroid energy and width of collective motion 
versus the pseudoparticle width in the low-temperature linear 
regime (|6^). The values are scaled to the free one a p = 0. 

In Figure |^ we give the ratio of the centroid energy and 
damping to the corresponding un-coarse grained ones for 
typical nuclear situation. We see that with increasing 
width of coarse graining the centroid energy decreases 
and the width increases. This is understandable because 
due to coarse graining we lower the particle- hole thresh- 
old resulting into lower centroid energy and an artificial 
damping at the same time. Since coarse graining is un- 
avoidable in numerical implementations the extraction of 
damping width of giant resonances should be critically 
revised with respect to the used pseudoparticle width. 

In contrast to the coarse graining in momentum space, 
we find a different behavior in spatial domain. In figure 
H we plot a realistic numerical solution of the full Vlasov 
equation describing monopole resonances in 40 Ca. Here 
the spatial width of the test particles has been varied. We 
see that the centroid energy is increasing with increasing 
testparticle width. This behavior can be understood from 
figure |]. The spatial coarse graining enlarges the gradi- 
ent of the potential and therefore increases the restoring 
force. This translates into a higher incompressibility and 
a higher collective energy. 
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FIG. 3. The centroid energy of monopole resonance in 40 Ca 
versus spatial testparticle width within a Vlasov simulation. 
The centroid energy is increasing with increasing width. 

A more complete discussion of the solution of Vlasov 
equation and also application to giant resonances far from 
the stability line can be found in pa] . 



mined nuclear potential is overestimated by testparticle 
simulations with finite width of the test particles. This 
should have implementations on Thomas- Fermi calcula- 
tions where an overbinding is found. 

The linear response function is calculated for homoge- 
neous systems and the spectra of density fluctuation is 
presented. It is found that the RPA polarization function 
becomes modified due to the finite momentum width of 
testparticles. These modifications can be understood as 
an internal structure the particles bear. This structure 
function is formally compared with vertex corrections to 
the RPA. It is pointed out that the higher order ver- 
tex corrections beyond RPA can be casted into similar 
structure functions. Therefore we suggest a method of 
simulating higher order correlations in one-body treat- 
ments by choosing an appropriate momentum width of 
testparticles. 

As a practical consequence the collective mode is an- 
alyzed. We find that the coarse graining enhances the 
damping width and lowers the centroid energy of col- 
lective modes, e.g. giant resonances. The dependence 
on the coarse graining width is given quantitatively and 
corresponding simulations should be revised. 
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VI. SUMMARY 

The dissipative features of coarse grained Vlasov equa- 
tions are investigated. This procedure occurs due to nu- 
merical simulation techniques. We have calculated ex- 
plicitly the entropy production, which is due to nonlin- 
ear mode coupling but not due to dissipation. We find 
that a sum of modified Boltzmann distributions is ap- 
proached by the coarse graining. Examples are shown 
where no stationary condition is approached at all since 
the solution oscillates and examples are given where a 
stationary solution can be reached. The different behav- 
ior is completely determined by the initial distribution 
and the used potential. The stability analysis leads to 
a criterion for the initial distribution dependent on the 
width of test particles, which can only lead to stationary 
solutions. 

We have demonstrated by a special model that the two 
steps, coarse graining and dynamical evolution of the dis- 
tribution function, are only interchangeable if the initial 
distribution is also coarse grained. It is argued that this 
property does not hold in general due to the feedback 
of the selfconsistent potential. Because the coarse grain- 
ing is unavoidable in numerically implementations, the 
Vlasov codes should be critically revised with respect to 
the question if they start really from a coarse grained 
initial distribution. 

Thermodynamical consequences are discussed. A cor- 
related part to any thermodynamical observable is calcu- 
lated explicitly. It is given in terms of the space and mo- 
mentum width. We find that the selfconsistently deter- 
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APPENDIX A: 



REFINED TESTPARTICLE 
PICTURE 



We want to consider now the question how the dy- 
namics is influenced by coarse graining. While the orig- 
inal Vlasov equation (|2|) is represented by the Hamilton 
dynamics of testparticles via (||), we like to know how 
the equation of motions are changed by coarse graining. 
Therefore we obtain with (||) the coarse grained kinetic 
energy K and the potential energy W as 



K(t) 



W(t) = - 



dpdr p 2 



(2tt^) 3 2m" 
dpdr p 2 + o-p 
(2tt^) 3 2m 
dpdr 
{2nh) 3 
dpdr 
{2nh) 



f{p,r,t) 
V s (r,t)f s (p,r,t) 
3 V{r, t)exp(a I 2 dr r )f(p, r, t). (Al) 



We see that the coarse graining in momentum space leads 
to an additional apparent temperature a 2 = mT, while 
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the coarse graining in space modifies the potential and 
therefore the dynamics. 

With the help of © and (g) we find from the varia- 
tion of the action J dt(K(t) — W(t)) just the Hamilton 
equations 



at m 



d 
dR 



■V s (R,t) 



(A2) 



irrespective which form we use of (Al). This is clear 
because the underlying dynamics is entirely covered by 
the Hamiltonian dynamics. Therefore even the center of 
mass movement of the Gaufi packets Ri , Pi follows the 
Hamilton dynamics. 

We may now turn the question around and define a 
new quasi(test)- particle picture. We seek the equation 
of motion for new testparticles (Pi, R) which should now 
represent the coarse grained distribution (^|) 

f(prt) = {fs} g 

OO 

= Y / 9r(r-R(t))g P (p-m)) 



2 = 1 
OO 



J2S(r-R(t))5(p-R(t)). 



(A3) 



»=i 



The last step is just the definition of the new test particles 
which obey the equation of motion derived from (Al) 



at m 



(A4) 



with V the coarse grained mean field (10). We observe 
that the coarse grained distribution / can be represented 
by a set of quasi(test)- particles which obey Hamilton 
equations with a modified effective potential 



V cB (r,t) =e-^ 9 *V(r,t) 

^V(r,t)-a 2 r d 2 r V(r,t). 



(A5) 



We see that the relation between the coarse grained mean 
field potential (0) and the effective one can be written 
as a double folding 



V(r,t) 



V eS (r,t) = {{V eS (r,t)} gr } gr 



V s (r,t)={V eS (r,t)} 6 



(A6) 



(A7) 



where we have used 



the relation {f(r)} gr — 
Please observe that the here 



exp(^^)/(r), see |6[ 
pre sented relation is just the inverse relation given in 
pq-Eo|. The difference comes from the inverse picture 



used here. The authors of |?8]-[l0| investigated how the 
dynamical equation change if in the action the distribu- 
tion functions are replaced by their coarse grained ones. 
We present here the opposite view that the action is un- 
changed but is rewritten according to coarse graining, 
l ead ing to modified explicit kinetic and potential energy 
(Al). This still leads to unmodified equation of motions 
for the center of mass coordinates of the Gaufi packets 
as outlined above. When we now represent the coarse 
grained distribution by a new set of sharp testparticles, 
the dynamics becomes modified and an effective potential 
appears which is the inverse folded mean field potential 

We conclude that we can map the coarse grained 
Vlasov equation to the original Vlasov equation if test- 
particles are introduced which obey a modified dynamical 
equation with the new inverse folded mean field poten- 
tial. Within the text we have followed the other route to 
describe the influence on the dynamics in the old picture 
to make the effects more transparent. 



APPENDIX B: COARSE GRAINING OF 
SELFCONSISTENT MODEL 



Here we give the explicit calculations of chapter III B 



a. First solving than coarse graining 

We like to choose as an initial condition a Gaussian 
distribution, which represents an equilibrium distribution 



/o(x) = (2n)- 3 A/Det{A}e-* 



with 



H 
A 



,-i x J Ax 



(Bl) 



(B2) 



Then the solution ( J33| ) of the selfconsistent Vlasov equa- 
tion reads 



/ s (x,t) = (27t)" 3 •v/Det{A} e ^^ Ax+c ' a < Ax+c 



= (2tt)- 3 ^Dct{A} e -^ xr + cr ( Ar r 1 ) ATAA ( x+A ~ lc ). 

(B3) 

We can now coarse grain the distribution with ( pl| ) to 
obtain the result 



/(x,t) = 



(27r) _3 e"^ Ax+c ) T ( A_1+:f:_1 ^ 1 ( Ax+c ' 



Det{(A T AA)- 1 +£} 



(B4) 



where 



12 



£ = ASA J 



(B5) 



with 



We see that the coarse grained solution (B4) can be 
represented at any time by the coarse grained initial dis- 
tribution 



M = (A T AA)" 1 - B. 



(B12) 



/(x,t) = {/ 5 (x,i)*S(x)}(y) 

= {/ (A X +C)*g(x,i)}(y) 
= {/o(x)*e(x,i)}(Ay + C) 



(B6) 



if we use the coarse graining ( |2l| ) but with the ( now time 
dependent) width (|B5|) . It is immediately obvious that 
the last identity in ( |B6[ ) is only valid if we use Gaussian 
initial distributions. Any other distributions will not al- 
low this rearrangement. Consequently the folding and 
the dynamics are not interchangea ble in general. 
The entropy is calculated from (B4) 



S = 31og(27re) + i logDet((A T AA)" 1 
with the explicit expression (|36|). 



(B7) 



b. First coarse graining than solving 

Equation (^|) is transformed into a partial differential 
equation of first order by Fourier transform x = (r, p) — > 
q = (q, x). The solution can be obtained for the Fourier 
transformed distribution / 

/(q,i) = /o((^ 1 ) T q)exp(iq T Bq + Dq) (B8) 

with the matrix A as defined in (p2f) and 



D = iva / dt'Q(t') 



sin(u;(t— t')) 
muj 

cos(w(i - t')) 



1 - CT. 



2 



B, 1— cos(2ujt) • , , . 

= —I — — — a VdUl) rntAj V / 

2 muj \ sin(2o>t) mu)(l — cos(2wi)) 



(B9) 

The selfconsistency requirement (^6|) leads to the same 
solution as ([mJ). This shows that the coarse graining 
does not affect the nonlinear feedback of the selfconsis- 
tent potential within this model. In other more nontrivial 
models this needs not to be the case. 

If we now use as initial distribution once more a Gaus- 



sian (Bl) which reads in Fourier transform 



/o(q) = cxp(— -q A l q) 



(BIO) 



we obtain the time dependent distribution after inverse 
Fourier transform 



/(x) 



(27T 



.-3 



'Det{M} 



: exp 



1 



(x-iD) r M- x (x-iD) 



(Bll) 
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